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Abstract. We solve the gauged Nambu-Jona-Lasinio model at leading 
order in the large N expansion by computing the anomalous dimensions of 
all the fields of the model and other gauge independent critical exponents by 
examining the scaling behaviour of the Schwinger Dyson equation. We then 
restrict to the three dimensional model and include a Chern Simons term 
to discover the ^-dependence of the same exponents where 9 is the Chern 
Simons coupling. 
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Recently there has been a resurgence of interest in examining the dy- 
namics of models with four fermi interactions due to the possible composite 
nature of the Higgs boson of the standard model, |l|, This is based partly 
on the observation that, similar to what occurs in the Gross Neveu model 
Hi, fermion pairs can bind together to form a bosonic particle whose mass 
is around twice that of one of the fermions. In [||, |2) it was discussed how 
this bound state can effect the spontaneous symmetry breaking of the Higgs 
boson in the more conventional understanding of the generation of mass 
for the particles we observe in the real world. To improve our knowledge 
of the basic four fermi model which underlies this alternative approach it 
is important to have a comprehensive picture of the quantum properties 
of such models. One technique widely used to probe these theories is the 
large N expansion. As a first step in this direction the O(N) Gross Neveu 
model has recently been solved at 0(1/N 2 ) in a 1/N approximation where 
N is the number of fundamental fermions, 0, H|. The method used was 
based on the earlier work of 0, |7| and involved computing the anomalous 
dimensions of the basic fields and bound states of the theory in arbitrary di- 
mensions by solving the theory near its d-dimensional critical point. There 
the fields obey a conformal symmetry and their critical exponents, which 
are related to the critical renormalization group functions, can be deduced 
by solving the Schwinger Dyson equations self consistently. More recently 
the chiral Gross Neveu or Nambu-Jona-Lasinio model, @, has also been 
solved at 0(1/N 2 ) in the same approach, [||. However, to fully understand 
the dynamics of these simplified models in relation to the standard model 
it is important to consider the more realistic scenario where four fermi in- 
teractions are coupled to gauge fields. In this letter we report the results of 
such an analysis for the Nambu-Jona-Lasinio model coupled to QED where 
we have a continuous global U(l) x U(l) chiral symmetry. We will deduce 
the anomalous dimensions of the fields in arbitrary dimensions in a large N 
expansion. This will be important for an e-expansion approach to under- 
standing which universality class the models lie in and for comparison with 
the perturbation theory of other models. 

As we will deal with a model with a U(l) gauge field it is also important 
to record new results for another area where four fermi models are also of 
interest. It has recently been argued that a model which underlies high 
T c superconductivity involves a four fermi interaction coupled to QED in 
three dimensions, |J. Whilst our results will produce three dimensional 
gauge independent exponents, which are therefore physically relevant, we 
will also couple our model explicitly to a Chern Simons term specifically in 
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three dimensions, [10-12]. The motivation for including such a topological 
term is its relation to anyons which are believed to be a mechanism for 
high T c superconductivity, |l^, |14|]. Therefore, it is important to produce 
gauge independent exponents which depend on the Chern Simons coupling 
9, which does not get renormalized in the quantum theory and therefore 



appears as a pure ^-independent parameter in the exponents, |15[. As the 
gauge independent exponents are physically interesting and measurable, this 
will allow one to determine whether 9 is non-zero or not and therefore if such 
topological terms exist in nature. 

The lagrangian of the theory we consider is 
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where a and ir are auxiliary bosonic fields whose elimination yields the usual 
four fermi interaction and which are fermion bound states in the quantum 
theory, §. We have scaled the coupling constant e and g in such a way 
that the couplings of the 3- vertices are unity and is our U (1) gauge field. 
The fields tp l lie in an iV-tuplet, 1 < i < N, and 1/N will be our expansion 
parameter. We follow the conventions of [16| in defining the properties of 



7 5 as {7^,7 } = and tr(7 5 7 w ...7^™) = 0, where we take trl = 4. 

To solve (1) via the self consistency approach of [4-7] the anomalous 
dimensions of the fields are deduced by substituting the asymptotic scaling 
forms of the fields at criticality into their skeleton Dyson equations with 
dressed propagators. More concretely for (1), we take in coordinate space, 
as x -> 0, [|, £7), 



ip(x) 



Aj 
(x 2 ) a 
D 



B 



a{x) ~ 



{x^y ' 

2p(l - b) 



7t(x) 



c 



2Vy 



(x 2 ) 



(2/x - 2p - 1 + b) 



(2) 



where b is the covariant gauge parameter. Retaining b in our calculation will 
prove useful in demonstrating which of our exponents are gauge independent 
and therefore physically relevant. The quantities A, B, C and D are the 
respective amplitudes of the fields and a, (3, 7 and p are the dimensions of 
the fields and we define their anomalous dimensions by 



o 



H + irj , P = 1 - ?7 - Xa , 7 = 1 - V ~ Xn , P = l-'H-XA (3) 
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where we define the spacetime dimension to be d = 2p, r] is the fermion 
anomalous dimension and x denotes the vertex anomalous dimension of the 
respective 3-vertices of (1) and they contain information on the structure 
of the quantum theory. To compute each of these, which depend on N, fx 
and are 0(1/N), one represents the skeleton Dyson equations with dressed 
propagators, which are valid in the critical region x — > 0, by the scaling 
forms (2). At leading order these are illustrated in fig. 1 and therefore they 
correspond to 

n , n , , , 2n[ /Q (l-6)-(/i-l)(2 / u-2p-l + b)] 

° = r(a-l) + * + y+ (2„ - 2p - 1 + ft) (4) 

= p(f3) + ANz (5) 
= p( 7 ) + ANy (6) 
(2p-2p-l)p(p) _ 4(a - l)u 

4N(n-p) (2a -1) { ) 

where z = A 2 B, y = A 2 C, u = A 2 D with r(a) = ap(a)/(p — a), p(a) 
= a(a — ^)/[7r 2/i a(a)] and a(a) = T(p — a)/T(a). Therefore eliminating 
z, y and u from (4)-(7) and substituting the leading order values for the 
exponents, one deduces 

r(2/i - 1) 



m 



r(i- M )r( M -i)/xr2( M ) 



(2/x-i)[(2 M -i)( At - 2) + M &] 1 



(8) 



4(/i - 1) 

where r/ = r ]il^ % an d (8) represents an analytic expresion which can be 
expanded in d = 2 + e or d = 4 — 2e to compare with other models. Whilst 
T)\ is 6-dependent it always must be computed first in any calculation since 
the other exponents are expressed as a function of it |J. In three dimensions 
we record 

4(36- 1) 

To determine the anomalous dimensions of the remaining fields one ex- 
amines instead of the 2-point function, the scaling behaviour of the respec- 
tive 3-point vertices in the critical region based on the method of jl8| adapted 
to the 4-fermi model with discrete chiral symmetry in Q. We find, however, 
that the extension of the Gross Neveu model with a continuous global chiral 
symmetry reveals a significantly different structure for the exponents com- 
pared to the coupling to the case of a discrete symmetry, [19]. For example, 
in arbitrary dimensions 

(2/x-l)(2 M -l + %r(2^-l) 

Xal XlTl " 4r( M + i)r 2 (/i)r(2- M ) 1 ' 
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xai = -m (ii) 

The absence of any 0(1/ N) correction to the gauge field exponent is merely 
a consequence of the U(l) Ward identity of the model which was observed in 
j pO[ , In the pure chiral Gross Neveu model, ||, Xai and Xwi are both zero 
in contrast to (10). The inclusion of the An field has given rise to quantum 
fluctuations which generate an anomalous piece at 0(1/N). Moreover, in 
(3 and 7 the gauge parameter b cancels in all dimensions so that we have 
deduced the first of the gauge independent exponents of this letter. 

The advantage of solving (1) via the self consistency approach of [4- 
7] is that it is relatively straightforward to push the analysis to the next 
order, 0(1/ N 2 ). This has been carried out separately for the 0(N) case, 
[||] , and for QED, where the fermion anomalous dimension has been 
computed at 0(1 /N 2 ) as well as in the amalgam model where the mixed 



graphs, which occur at next order, were computed [22]. We have been able 
to extend that analysis to (1). This involves including the higher order 
graphs in the Dyson equations and we have illustrated a representative set 
for each Dyson equation in fig. 2. There are three higher order two loop 
graphs for each of the cr, n and equations, whilst there are nine for ip. 
Due to the high degree of symmetry present the values of only two of these 
are independent. Restricting to the Landau gauge we determine, 

_ r ( y - do, - if . ( „ - - 1) /» 

12 "[ 2 4 \(P-1) (M 

+ 3 ^~_ 1) 1) [(2f " 1) " (f " I) 2 + 2(« - 1) 4 (2(. + 1)(2^ - 3)] 

+ l{2 " -Dfr-5 + 'b-i)1 ( 6/ / ( „ _ im _ i) fe - 1 



32fx V V (/J. — 1) 

(12 



_ 2^-2)(2,-1) +(2m _ i)2(m _ 2)+4(m _ i)2 

where ^ = T(2fi- 1)/[T(2 - fi)T 2 (fi)T (fi + 1)], 6(/i) =^'(//-l) - ^'(1) and 
^(n) = ip(2fj, — 2) + ip(3 — //) — -0(1) — ^(m) and ip(x) is the logarithmic 
derivative of the T-function. The result (12) represents the first 0(1/N 2 ) 
expression available for (1) and in three dimensions, 

4[27vr 2 + 14] 

m = 27 ,4 (13) 

To complete the 0(l/iV) analysis of (1) we need to compute the correc- 
tions to the exponent A = /x — 1 + 0(l/iV), which has been studied in other 



5 



models, [Q, 17, [HJ. The first step in the method to achieve this involves 
considering corrections to the asymptotic scaling functions (2) of the form 



vr(x) ~ ^_[l + C'(x 2 ) A ] (14) 



D 



2p(l - 6) 



(2/m - 2p - 1 + 6) x 2 
, D , a / „ , 2(p-A)(l-6) x M x, 



within the critical Dyson equations, where A', B' , C and D' are new ampli- 
tudes. This procedure is more laborious than for the original formulation of 
[|J since, as was noted in llTj, one has to include corrections coming from the 
a, it and two loop self energy graphs, where there is an (x 2 )^ insertion 
on the respective lines. Omitting these pieces would result in an erroneous 
analysis. To determine the 0(1/N) correction to A the consistency equa- 
tion based on (14) decouple into a set for 772 already considered and a set 
involving A' , B 1 , C and D'. To ensure this is a consistent set of equations, 
the determinant of the matrix formed by these amplitudes as basis vectors 
must vanish. With the explicit values of the 2-loop integrals calculated in 



[ 19 1 , the determinant yields two solutions. We discard the trivial one Ai = 
0, in favour of 

Ax = - ±(2/x-1)(2/j 2 - 2/i + 1)7/1 (15) 

and we have checked that the 6-depedendence has cancelled. 

We finish the analysis by restricting our arbitrary dimensional work to 



three dimensions and include the term, [11, 12 1, 



Oe^A^A" (16) 

in (1), where we work in euclidean space, £123 = 1 and 9 is the Chern Simons 
coupling constant. (The (F^ v ) 2 term of (1) is now omitted and the transverse 
contribution to the A^ propagator will be generated in the large N expansion 
through the relevant graphs of fig. 1.) The aim will now be to deduce the 
^-dependence of the exponents calculated at 0(1/N) due to the potential 
relation of four fermi models to high T c superconductivity. As was noted in 
[|i"9f one advantage of calculating with (16) in the large N self consistency 
approach of || [7j is that the spacetime dimension is fixed and therefore 
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one does not need to extend the definition of e^ ua to d-dimensions as one 



would in a perturbative analysis using dimensional regularization, [23]. The 
only consequence that (16) has on the scaling forms is the appearance of 
an overall factor (1 + # 2 ) -1 in the scaling form of the gauge field two point 
function A~^,(x), [jOJ. With this observation it is straightforward to repeat 
the 0(1/N) analysis for 6 ^ in d = 3 to find 

4[36 - 1 + 8 2 } 4(6 + 2) 



3tt 2 (i + 2 ) ' A<T1 ™ x n 2 (i + e 2 ) 

8{9 2 + 5) 

XAi = - VI , Ax = - 3?r2(1 + g2) (17) 

where each exponent is a monotonic function of 8. Similar 8 dependent 
exponents have been obtained at leading order in 1/N in other models, 



[P0| , 24, 25 1 . As a check on the correctness of (17), they agree with the 
earlier 8 = values, (9), (10) and (15), and as 8 — > oo one recovers the 
previous results of || since then the gauge field decouples from (1) and (16) 
where clearly one is left simply with the pure Nambu-Jona-Lasinio model. 
This latter limit has also been examined in other models, such as the bosonic 
a model on CP{N) [^4], pEfl , and the same decoupling feature of the gauge 
field is observed there as well. 

We conclude with several remarks. First, we have solved the gauged 
Nambu-Jona-Lasinio model or chiral Gross Neveu model at leading order 
large N by deducing the gauge independent exponents of the model. As 
these expressions are in arbitrary dimensions they will be of use in deter- 
mining agreement with the e-expansion of the explicit perturbative renor- 
malization group functions of the theory. By including a Chern Simons 
term to reveal non-trivial 8 dependence, these results will be useful in the 
experimental situation of differentiating between the nature of the four fermi 
interaction since the 8 dependence of our results, (17), differ substantially 



from those of IS]. Finally, it is worth pointing out that the conformal 
methods we have used here to deduce exponents are much more efficient in 
terms of computing integrals than the conventional large iV methods and, 
moreover, have the advantage of allowing one to easily probe the quantum 
structure to next to leading order. 



7 



References. 

[1] W.A. Bardeen, C.T. Hill & M. Lindner, Phys. Rev. D41 (1990), 1647. 

[2] A. Hasenfratz, P. Hasenfratz, K. Jansen, J. Kuti & Y. Shen, Nucl. Phys. 
B365 (1991), 79. 

[3] D. Gross & A. Neveu, Phys. Rev. D10 (1974), 3235; Y. Nambu & G. 
Jona-Lasinio, Phys. Rev. 122 (1961), 345. 

[4] J.A. Gracey, Int. J. Mod. Phys. A7 (1991), 395; 2755(E). 

[5] J.A. Gracey, Phys. Lett. 297B (1992), 293. 

[6] A.N. Vasil'ev, Yu.M. Pis'mak & J.R. Honkonen, Theor. Math. Phys. 

46 (1981), 157. 

[7] A.N. Vasil'ev, Yu.M. Pis'mak & J.R. Honkonen, Theor. Math. Phys. 

47 (1981), 291. 

[8] J.A. Gracey, 'The Nambu- Jona-Lasinio model at 0(l/iV 2 )', LTH-308 
to appear in Phys. Lett. B. 

[9] A. Kovner & D. Eliezer, Int. J. Mod. Phys. A7 (1992), 2755. 

[10] R. Jackiw & S. Templeton, Phys. Rev. D23 (1981), 2291. 

[11] J. Schonfeld, Nucl. Phys. B185 (1981), 157. 

[12] S. Deser, R. Jackiw & S. Templeton, Phys. Rev. Lett. 48 (1982), 975. 

[13] E. Witten, Commun. Math. Phys. 121 (1989), 351. 

[14] I. Dzyaloshinsky, A.M. Polyakov & P. Wiegman, Phys. Lett. 121A 
(1987), 112. 

[15] A. Blasi, N. Maggiore & S.P. Sorella, Phys. Lett. 285B (1992), 54. 

[16] S. Hands, A. Kocic & J.B. Kogut, Ann. Phys. 224 (1993), 29. 

[17] J.A. Gracey, 'Algorithm for computing the /3-function of quantum elec- 
trodynamics in the large Nf expansion', LTH-294 to appear in Int. J. 
Mod. Phys. A. 



8 



[18] A.N. Vasil'ev & M.Yu. Nalimov, Theor. Math. Phys. 55 (1983), 423; 
56 (1983), 643. 

[19] J.A. Gracey, Ann. Phys. 224 (1993), 275. 

[20] J.A. Gracey, J. Phys. A25 (1992), L109. 

[21] J.A. Gracey, Mod. Phys. Lett. A7 (1992), 1945. 

[22] J.A. Gracey, J. Phys. A26 (1993), 1431. 

[23] W. Chen, G.W. Semenoff & Y.S. Wu, Phys. Rev. D44 (1991), R1625; 
L.V. Avdeev, G.V. Grigorev & D.I. Kazakov, Nucl. Phys. B382 (1992), 
561. 

[24] S.H. Park, Mod. Phys. Lett. A7 (1992), 1579. 

[25] G. Ferretti & S.G. Rajeev, Mod. Phys. Lett. A7 (1992), 2087. 



9 



Figure Captions. 

Fig. 1. Leading order skeleton Dyson equations. 

Fig. 2. Basic structure of graphs for 0(1/N 2 ) corrections. 
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